Extragalactic Jets: Theory and Observation from Radio to Gamma Ray 

ASP Conference Series, Vol. , 2007 

T. A. Rector and D. S. De Young (eds.) 



Electron Energy Distributions at Relativistic Shock Sites: 
Observational Constraints from the Cygnus A Hotspots 



C. C. Cheung 1 , L. Stawarz 

Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University, Stanford, CA 94305, USA 

D. E Harris 

Harvard-Smithsonian Center for Astrophysics, 
60 Garden St., Cambridge, MA 02138, USA 

M. Ostrowski 

Astronomical Observatory, Jagiellonian University, 
ul. Orla 171, 30-244 Krakow, Poland 

Abstract. We report new detections of the hotspots in Cygnus A at 4.5 
and 8.0 microns with the Spitzer Space Telescope. Together with detailed pub- 
lished radio observations and synchrotron self-Compton modeling of previous 
X-ray detections, we reconstruct the underlying electron energy spectra of the 
two brightest hotspots (A and D). The low-energy portion of the electron dis- 
tributions have flat power-law slopes (s ~ 1.5) up to the break energy which 
corresponds almost exactly to the mass ratio between protons and electrons; we 
argue that these features are most likely intrinsic rather than due to absorp- 
tion effects. Beyond the break, the electron spectra continue to higher energies 
with very steep slopes s >3. Thus, there is no evidence for the 'canonical' s=2 
slope expected in 1st order Fermi-type shocks within the whole observable elec- 
tron energy range. We discuss the significance of these observations and the 
insight offered into high-energy particle acceleration processes in mildly rela- 
tivistic shocks. 



1. Introduction 

Cygnus A is the nearest example (z=0.056) of a powerful Fanaroff-Riley type- 
II radio galaxy. As such, it contains some of the brightest radio hotspots to 
study these classical 'working surfaces' of relativistic jets (Hargrave & Ryle 1974; 
Blandford & Rees 1974). Although the hotspots are well-studied at radio and 
X-ray wavelengths (§ 2), Cygnus A is viewed through the galactic plane, making 
them difficult to study at optical wavelengths (1.26 mag of extinction at V-band; 
Meisenheimer et al. 1997). Additionally, it is well known that bright foreground 
stars lie in the direction of Cygnus A, and this has also hindered optical studies 
of the hotspots. To circumvent these obstacles, we obtained new Spitzer Space 
Telescope images (Figure 1) resulting in successful mid-infrared detections of 
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the hotspots. These new data, together with a wealth of published high-quality 
measurements at other frequencies, give us some of the most complete hotspot 
SEDs for a powerful radio galaxy. This contribution is a summary of our recently 
published study (Stawarz et al. 2007). 



Cygnus A: Spitzer & VLA 





Figure 1. Spitzer (our work) and VLA (green; from Perley et al. 1984) 
image of Cygnus A. The 4.5/mi sources are shown in blue, 8.0/im emission in 
red, and those emitting in both IR bands in pink. 



2. Spitzer Imaging and the Hotspot SEDs 



We obtained new Spitzer /IRAC ( Fazio et al. 20041 ) observations simultaneously 
at 4.5 and 8.0 //m for 3.6 hrs. Galactic extinction amounts to only 0.02-0.03 
mag at these bands. A foreground star is known to lie within an arcsecond of 
the position of hotspot A - extrapolating its spectrum into the IRAC bands, its 
contamination to the hotspot IR fluxes is expected to be negligible (see Fig. 2g 
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in Meisenheimer et al. 1997). Radio data were obtaine d from published VLA 



and BIMA measureme nts and maps ( Carilli et al. 1991 ; Wright &: Birkinshawl 
120041 : lLazio et allfeOOd ). We focus our discussion on the two brightest hotspots 
(A and D), whose spectral energy distributions (SEDs) are shown in Figure [2j 




Figure 2. SEDs of the two brightest radio hotspots (A, in the western lobe 
and D in the eastern one) in Cygnus A. The data points are shown with error 
bars and open triangles are upper limits. The solid lines are our fits to the 
synchrotron and SSC components of the SEDs. Dashed lines show the effect 
of adding an additional break in the synchrotron spectrum (i^, r = 0.5 x 10 12 Hz 
for D and 1.2xl0 12 Hz for A). The dotted lines correspond to a hypothetical 
increase of a\ to 0.5. 

We modeled the synchrotron (radio-to-infrared) portion of the hotspot SEDs 
with a double power-law that decays exponentially at the higher frequencies: 

gsyn ^ f v~ ai for i/ min < v < v ct ,^ 

v \ u- a2 exp (-v/v max ) for v > u CI 

The fitted parameters (cti, «2, ^cr, ^max) are (0.28, 1.2, 2.6xl0 9 Hz, >3.3xl0 13 
Hz) for hotspot A, and (0.21, 1.1, 3.1 xlO 9 Hz, 0.9xl0 13 Hz) for hotspot D. The 
minimum frequency, z/ m j n , can be extrapolated down to 74 MHz with the VLA 
detections of the hotspots by Lazio et al. (2006) at this frequency. 

The hotspots emit strong X-ray emission by the synchrotron self-Compton 
(SSC) process as discovered with ROSAT (Harris et al. 1994). Assuming this for 
the Chandra detection of the hotspots (Wilson et al. 2000; Wright & Birkinshaw 
2004), we infer magnetic field strengths of ~0.2-0.3 mG. The emitting plasma is 
close to equipartition with U e /UB = 3—8 (Figure 3), consistent with the previous 
fits of the Chandra data. An HST optical measurement of hotspot D (Nilsson et 
al. 1997) lies near the low energy extrapolation of the SSC spectrum; this makes 
it one of the few hotspots with detected optical SSC emission (Brunetti 2004). 
Increasing the low energy spectral index {a\j of hotspot D from its observed 
value to the 'canonical' value of 0.5 would overproduce the observed optical and 
IR emissions, so the observed flat spectra are likely intrinsic. 

With the magnetic field and observed spectral fits, we can reconstruct the 
underlying electron energy spectra, n e (7) oc (s = 2a + 1; Figure 3). They 
have the following features: (i) The critical break energy corresponds almost 
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Figure 3. [left] Ratio of energy densities in ultrarelativistic electrons and 
magnetic field for the two brightest radio hotspots as a function of the low 
energy cutoff in n G (7). Solid lines correspond to the case with observed low- 
frequency synchrotron spectral indices of ot\ = 0.21 (hotspot D; lower line) and 
0.28 (hotspot A; upper line), while dotted lines correspond to a hypothetical 
increase of ol\ to 0.5 for both features, [right] Electron energy spectrum 
for hotspots A (dashed line) and D (solid line), as inferred from the SSC 
modeling. Shaded regions show the energy ranges affected by absorption or 
radiative losses. 



exactly to the ratio of the proton to electron masses, 7 cr ~ m p /m e . (ii) The 
low-energy segment of the electron distributions continues down to at least ~0.1 
m p /m e , with flat spectra; the power-law slopes are s ~l.Jf.-L6. (Hi) The new IR 
data constrain the high-energy electron spectra, with indices s>3 and maximum 
energies ^ 50 m p /m e . The "standard" electron spectrum expected from the 
diffusive (1st order Fermi) shock acceleration in the non-relativistic test-particle 
limit, n c {j) oc 7~ 2 , is not observed. The overall curved spectra are likely to be 
intrinsic rather than due to any absorption/cooling effects. 



3. The Cygnus A Hotspots as Relativistic Shocks 

Terminal shocks of powerful jets, like those in Cygnus A, are probably (mildly) 
relativistic, with an oblique magnetic field configuration. As such, it should 
not be surprising that the diffusive acceleration process known from the non- 
relativistic test-particle models are not directly applicable, i.e., the lack of a 
canonical shock- type spectrum oc 7~ 2 - see Begelman & Kirk (1990) and Niemiec 
et al. (2006) for relevant discussions. The break frequency corresponding to the 
m p /m e mass ratio preceded by a fiat-spectrum power-law is also an interest- 
ing finding, which confirms some previous evidences for such in other powerful 
hotspots (see Leahy et al. 1989). 

Our results, although inconsistent with the simplest modeling of terminal 
hotspots in radio galaxies as sites of the efficient 1st order Fermi acceleration, are 
instead in good agreement with our current understanding of relativistic shock 
waves. In particular, a flat power-law electron spectrum up to m p /m e energies is 
consistent with the simple (ID) resonant acceleration model discussed by Amato 
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& Arons (2006) for relativistic shocks. Their simulations show that as a result 
of resonant absorption of cyclotron emission produced by cold protons reflected 
from the shock front, a flat power-law electron energy distribution emerges be- 
tween energies 7 ~ T and T(m p /m e ), the latter possibly reduced by a factor 
of a few due to thermal dispersion in the upstream proton momenta (T is the 
bulk Lorentz factor of the upstream medium). The slope of this distribution is 
not well constrained by the simulations, but can be relatively flat and appears 
to depend on the plasma content (i.e., on the ratio of proton to electron num- 
ber density). Although these particular simulations involve a large T plasma, 
analogous processes are expected to take place also at mildly relativistic shocks 
considered in hotspots like Cygnus A (r ~ few), so the model results can at 
least be applied qualitatively. 

If such an association is correct, and if the results of Amato & Arons (2006) 
are applicable, then it would automatically imply that (i) the jets in powerful 
radio sources like Cygnus A are made of an electron-proton rather than electron- 
positron plasma, (ii) a significant fraction of the jet kinetic power is carried by 
cold protons, and (Hi) the number density of protons within the jets is most likely 
lower than that of the electron-positron pairs (see the discussion in Stawarz et 
al. 2007). We note that these conclusions regarding the jet content would then 
be consistent with those presented for quasar jets (e.g., Sikora & Madejski 2000), 
which are widely regarded as beamed counterparts of FR II radio galaxies like 
Cygnus A. 
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